Dissipative dynamics of circuit-QED in the mesoscopic regime 
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We investigate the behavior of a circuit QED device when the resonator is initially populated 
with a mesoscopic coherent field. The strong coupling between the cavity and the qubit produces 
an entangled state involving mesoscopic quasi-pointer states with respect to cavity dissipation. The 
overlap of the associated field components results in collapse and revivals for the Rabi oscillation. 
Although qubit relaxation and dephasing do not preserve these states, a simple analytical description 
of the dissipative dynamics of the circuit QED device including cavity relaxation as well as qubit 
dissipation is obtained from the Monte-Carlo approach. Explicit predictions for the spontaneous 
and induced Rabi oscillation signals are derived and sucessfully compared with exact calculations. 
We show that these interesting effects could be observed with a 10 photon field in forthcoming 
circuit QED experiments. 

PACS numbers: 03.67.Mn,03.65.Yz,32.80.-t,74.50.+r 



I. INTRODUCTION 

Cavity quantum electrodynamics is an ideal play- 
ground for the realization of tests of quantum processes, 
simple quantum information processing and exploration 
of the quantum/classical boundary. Mesoscopic superpo- 
sitions made of coherent field components with different 
classical attribucs (phase and amplitude) and their decay 
as a result of decoherence have been observed in exper- 
iments involing circular Rydberg atoms [l[. Recently, 
an experimental scheme for cavity QED experiments us- 
ing Joscpshon qubits embedded in a superconducting mi- 
crostrip planar resonator has been proposed Q and real- 
ized Q by the Schoclkopf group at Yale. Relying on con- 
ventional lithography and nanofrabrication techniques, 
these circuit-QED devices open the way to scalable quan- 
tum circuits coupling qubits to high-Q cavities acting as 
quantum buses. Recent experimental progresses on these 
devices have been dramatic, leading to the experimental 
demonstration of quantum non-demolition measurements 
of the qubit state in the dispersive regime [J and of the 
resolution of photon number states of the cavity @ . 

These recent developements provide a strong motiva- 
tion toward studying the dynamics of circuit QED de- 
vices at the quantum/classical boundary. In particular, 
it is important to analyze their ability to produce meso- 
scopic Schrodinger cat states involving coherent compo- 
nents with a mesoscopic number of photons. Such states 
can be produced in the dispersive regime where the off 
resonant qubit acts as a state dependant transparent di- 
electrics. In this regime, the interaction between a coher- 
ent field and a qubit, initially prepared in a superposition 
of states naturally leads to a quantum superposition of 
phase shifts. 

A much faster way to produce entangled qubit + cavity 
states is to use a resonant qubit + cavity interaction. In 
the mesoscopic regime, photon graininess rapidely casts 
the initially coherent cavity field into a superposition of 
two components with different phases [3]. This phase 



splitting is a mesoscopic effect that disappears in the 
classical limit of a very large field, which is then left un- 
affected by the atoms. Thus, in the mesoscopic limit, the 
cavity field acts as a which-path detector for the atomic 
states. The overlap of field components of the qubit + 
cavity state is responsible for the collapse and revivals of 
Rabi oscillations which can be either spontaneous or in- 
duced by an echo sequence. This resonant phase splitting- 
effect has been evidenced in Rydberg atom experiments 
for fields containing up to a few tens of photons Its 
coherence has been checked using an echo technique bor- 
rowed from NMR Jj|, following a proposal by Morigi et 
al 0. 

In the presence of dissipation, the entangled qubit + 
cavity state is subject to decoherence which leads to a 
damping of the spontaneous and induced Rabi oscillation 
revivals. In a previous work [13], we have developped a 
simple analytical model that describes the behaviour of 
N identical non dissipative qubits resonantly and sym- 
metrically coupled to a high-Q cavity. Considering cav- 
ity relaxation as the source for decoherence, this model is 
perfectly well suited for describing Rydberg atom exper- 
iments. But this is not sufficient for circuit QED since it 
does not take into account qubit relaxation and dephas- 
ing. 

The purpose of the present paper is to present a simple 
analytical model that takes into account qubit relaxation 
and dephasing for a single qubit in a cavity thus general- 
izing results previously obtained by Gea-Banacloche 
and ourselves [Io| . Using the physical insight provided by 
the stochastic wave function approach [12j |. an analytic 
description for the decoherence of the mesoscopic qubit + 
cavity state is derived. Our derivation also sheds light on 
the range of validity of our previous analysis: whenever 
the internal interactions within a mesoscopic device tend 
to produce pointer states with respect to the dominant 
coupling to the environment, its evolution can described 
in term of these pointer states and decoherence function- 
als reflecting the cumulative imprints they left in the en- 
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vironment. In circuit-QED devices, this simple image is 
broken by qubit dissipative processes but a careful analy- 
sis of the dynamics enables us to obtain simple analytical 
results for the spontaneous and induced Rabi oscillation 
signals. In principle our analysis can be also used to dis- 
cuss field tomography and can also be extended to the 
case of several qubits. 

This paper is organized as follows. In Sec. HH the ba- 
sic model for cavity-QED is presented, its dynamics in 
both resonant and dispersive regimes are briefly recalled. 
In Sec. IIIIl circuit-QED devices are presented and dis- 
sipation sources are discussed. In Sec. IIV1 the dissipa- 
tive dynamics of a one qubit circuit-QED device is stud- 
ied within the stochastic wave function framework. For 
completeness, analytical results are derived both for the 
resonant and dispersive regimes, recovering in the latter 
case results from the Yale group Section fVl presents 
numerical results obtained from quantum Monte-Carlo 
simulations in the resonant regime. These are used to 
discuss the validity of our analytical model and to derive 
experimentally accessible windows for the observation of 
mesoscopic qubit + cavity states in forthcoming circuit- 
QED devices. In the conclusion, we comment on the 
inclusion of other dissipative effects which may be rele- 
vant in the study of other cavity QED superconducting 
circuits 



cigenstate are very close to the ones of the uncoupled 
cavity + qubit states. Performing a second order expan- 
sion in g/A leads to 

ft _1 -Hdis P = (wo + XC Z ) a) a + i(w qb + x) <** ■ (2) 

where x = 3 2 /4A represent the ac-Stark shift per pho- 
ton. In the presence of a mesoscopic coherent state with 
average photon number n, the cavity + qubit coupling 
is enhanced by the coherent field. Off diagonal terms in 
the Jaynes Cummings multiplet are small compared to 
diagonal ones when |A| 3> g\ffi. 

In the dispersive regime, qubit flips can only occur be- 
cause of qubit relaxation or external driving of the sys- 
tem. With respect to the cavity mode, the qubit behaves 
as a transparent medium with a state-dependent refrac- 
tion index. In Rydbcrg atom experiments, this regime 
has been used to produce Schrodinger cat states Ell, to 
perform QND measurement of the photon number [la ] 
and a measurement of the Wigner function of the field 
[T^ . In circuit QED devices, the cavity has been used 
to performed a QND measurement of the qubit state 
More recently, the strong coupling dispersive regime 
where a single photon drastically alters the qubit absorb- 
tion spectrum has been has been studied theoretically 
fl3| and demonstrated experimentally [j| . 



II. CAVITY QED 

A. The model 

The effective Hamiltonian describing the cavity QED 
systems is the Jaynes- Cummings Hamiltonian [l5| in- 
volving the coupling of an harmonic mode to a two level 
system (the qubit): 



/T^jc = loo a f a + ^ ° z + ?(oV~ + a a+) 



(1) 



In the strong coupling regime, the coupling energy g 
is assumed to be much larger than all energy scales 
characterizing dissipative processes in the system, i.e. 
the cavity relaxation rate re and the atomic relaxation 
and dephasing rates. Neglecting damping, the Jayncs- 
Cummings Hamiltonian can be exactly diagonalizcd since 
the Hilbert space decomposes into two dimensional mul- 
tiplets generated by the states |e, n) = |+) ® \n) and 
\g,n+l) = | — ) <g> \n). As a function of the qubit / cavity 
detuning A = w q b — i^o> the cavity + qubit system has 
two different regimes. 



B. The dispersive regime 

The dispersive regime is reached when the cavity and 
the qubit are out of resonance. For an initally empty 
cavity, it is obtained for |A| 3> g. In this regime, the 



C. The resonant regime 

The resonant regime is obtained when A = 0. In 
this case, the Jaynes-Cummings eigenstates are symmet- 
ric and antisymmetric combinations of the form (|e, n) ± 
+ l))/v / 2 separared by an energy g (vacuum Rabi 
splitting). This leads to vacuum Rabi oscillations be- 
tween states |e,0) and \g, 1) which have been observed 
18 1 and used to transfer the qubit state in the cavity 
19 ], In the mesoscopic regime n 3> 1, the resonant in- 
teraction leads to an entangled atom-field state with two 
quasi-coherent field components with different classical 
phases @ . The approximate solution by Gea-Banacloche 
precisely provides explicit expression for these states and 
provides the basic framework for discussing the complete 
dissipative dynamics of circuit QED devices. 

It has proven very convenient to describe the dy- 
namics using an effective Hamiltonian approach that in- 
volves effective spin operators acting within the Jaynes- 
Cummings multiplet fill [20| . Within the framework of 
this mesoscopic approximation (n 3> 1 and t <C n/g), 
the evolution of the cQED system can be computed ex- 
actly. Under the Jaynes-Cummings Hamiltonian, in the 
mesoscopic regime, the states 
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where \X^>) = (\e,p)±\g,p+l))/V2 remain of the same 
form with a time dependent angle 8(t) = g. In the meso- 
scopic approximation, \fy±(9)) can be approximated by 
a factorized state of the form [l(| : 

|**(0)) ~e^^ 2 |I} ± (0))®|^±(0)) (4) 
where the atomic polarization is given by: 

|D ± (0)) = -l(± e ^| + } + |-» (5) 
and the field component is 

|V±(0)> = e ±ie ^ 2 e-"/ 2 eT ^ /2 \ k ) • ( 6 ) 

The field state (0)) is proportional to a quasi coherent 
state whose parameter is roughly equal to \ffi e Tlt ^ 4v/ ^ *. 
At fixed n, the phase, in Fresnel plane, associated with 
the generalized Gea-Banacloche state is </>,„, e = 

—mQj2\ffi. The atomic polarization evolves in the equa- 
torial plane of the Bloch spere, its phase being perfectly 
correlated to the electromagnetic phase. 

The states \ip± (8)) act as a path detector for the atomic 
polarizations \D±{9)) and this explains the pattern of col- 
lapses and revivals of Rabi oscillations. When t -C g~ l , 
the field components \if>±(gt)} still overlap and therefore, 
Rabi oscillations arising from interference between the 
\D + ) and \D-) are observed. Once gt > 2ir, the field 
components do not overlap and Rabi oscillations disap- 
pear. Only when gtj1\ffi is a multiple of 27r, the two 
states overlap again. The information about the qubit 
state stored by into the cavity is erased, thus leading to 
a revival of Rabi oscillations. Henceforth, collapses and 
revivals of Rabi oscillations in the mesoscopic regime are 
a direct illustration of the complementarity principle. 

When gt/2^/E = ir, the qubit polarizations will coin- 
cide. Then, the atom + cavity state disentangles and 
the resulting cavity state is a mesoscopic Schrodinger cat 
states involving two quasi-coherent components with op- 
posite phases. This occurs much faster than in the dis- 
persive regime |A| 3> gVfi where such a Schrodinger cat 
state would be generated in 7r|A|/g 2 ^> 2iry/E/g. 



III. CIRCUIT QED DEVICES 

Recent experiments performed at Yale involve the cou- 
pling of a superconducting Joscphson qubit and a super- 
conducting microstrip resonator with high quality factor 
Q ■ We will first present these devices and recall how they 
can be described by the Jaynes-Cummings Hamiltonian. 
Dissipation sources will then be discussed. 



A. The cavity: a superconducting microstrip 
resonator 

The resonator is ID a transmission line made by 
photolithography whose lowest mode lies within the 1- 
10 GHz frequency range and whose quality factors can 
be as high as 10 6 [2l|. Besides this, the electromagnetic 
field associated with their eigenmodes is confined within 
a relatively small volume, thus leading to high electric 
fields between the center and ground planes (typically 
0.2 V/m). Such a resonator is usually characterized by a 
wave velocity v, a real impedance R ans its length L. 

The finite length superconducting resonator has many 
stationary modes up to a certain high frequency cutoff 
but in the present situation, only the coupling to the 
lowest energy mode The voltage difference between the 
inner and outer electrodes at point x is then expressed in 
terms of the mode creation and destruction operators: 

V{x)=V Q (p + ct) sm(^) (7) 

where Vo denotes the maximum voltage felt by the qubit. 

B. The artificial atom: a Josephson qubit 

Within the resonator lies a superconducting Joseph- 
son charge qubit (initially a Cooper pair box) coupled 
to the cavity electric field. It consists in a small super- 
conducting island connected to a superconducting reser- 
voir through a very thin insulating barrier. When the 
temperature ksT and the the Coulomb charging energy 
of the island E c are smaller than the superconducting 
gap, the qubit degrees of freedom are encoded by the 
charge state of the superconducting island. This de- 
vice can be characterized by the Josephson amplitude 
Ej and the Coulomb energy E c . Control parameters arc 
the voltage potential V g imposed to the superconducting 
island and the Josephson amplitude which can be tuned 
using a SQUID instead of a single Josephson junction 
[22l ]: Ej = Ej cos (27T$j/$ ) where $ = 2e/h denotes 
the flux quantum and <!>j the magnetic flux through the 
SQUID. 

In the charge regime E c <C Ej and for gate charge 
rig = C g Vg/2e restricted to a unit charge interval < 
n g < 1, the system effectively behaves like a two level 
system (TLS) involving two adjacent charge states of the 
island. In this charge basis, its effective Hamiltonian is 
given by: 

ffqb = -^-^0* (8) 

where B z = 4E C (1 — 2n g ) and B x = Ej. 

C. The qubit /cavity coupling 

The electromagnetic mode of the ID cavity naturally 
couples to the qubit charge through capacitance between 
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the resonator inner and outer electrodes. The coupling g 
energy is limited by the dipolar energy of the qubit within 
the resonator: g/uio < X c y/R/R q where A c is a capaci- 
tance ratio and R q = h/e 2 is the quantum of resistance. 
Classical electrodynamics shows that R is of the order of 
2R q ak/ e r where a is the fine structure constant, e r is 
the relative permittivity and A a geometric factor, usu- 
ally logarithmic in the microstrip aspect ratios. Within 
the 1 to 10 GHz frequency range, values of g/2n range 
from 11 to 105 MHz (see table [B the latter value being 
obtained using the "transmon" [23]. Assuming that the 
two level description of the qubit is valid, the coupling 
Hamiltonian between the qubit and the cavity is given 
by: 

H int = hg(a + at)(l-2n g -a x ). (9) 

The coupling being small compared to the resonator 
and qubit eigenfrcqucncies, the rotating wave approxi- 
mation is valid and close to the charge degeneracy point 
n g ~ 1/2, cq. ^ reduces to the Jaynes-Cummings 
Hamiltonian (fT|). 

D. Measurement protocols 

In these experiments, the Josephson qubit cannot be 
measured directly. This circuit-QED system is probed 
through the resonator which is then connected to external 
ports at its ends. A network analyzer is used to analyzed 
the phase and amplitude of a classical electromagnetic 
wave transmitted through the device. This method has 
been used to perform a non-destrictive measurement of 
the qubit state using a dispersive measurement technique 
Q. Thus, the relaxation time scale of the resonator kT 1 
sets a lower bound on the measurement time. Probing 
the state of the qubit thus requires that k > 72 where 
72 = 71 + 7 y /2. Although cavities with quality factors 
of the order of 10 6 have been manufactured [2l|, devices 
used in experiments have a lower Q in order to satisfy 
the fast measurement constraint. 

In order to avoid this limitation, one has to rely on an 
alternative method for measuring the qubit state. Re- 
cently, a new detection scheme based on the dynami- 
cal bifurcation of Josephson junction has been realized 
and provides a high contrast, low backaction, dissipation- 
less and fast measurement of the qubit state [24|. Thus, 
rapid improvement of Joscpshon circuit technology sug- 
gests that alternative methods of detection might become 
available in the near future. 

In the present paper, we focus on the evolution of the 
qubit + cavity system prior to measurement. In partic- 
ular, our main objective is to describe the dynamics of 
the circuit QED system at resonance in the mesoscopic 
regime and shed light on the main physical effects in- 
depcndantly from the measurement protocol. Therefore, 
for the sake of simplicity, we will not attempt modeling 
the subsequent measurement process. 



E. Dissipation mechanisms 

1. Qubit relaxation and dephasing 

Qubit dissipation and decoherence arise from their 
coupling to environmental degrees of freedom either ex- 
trinsinc (measurement circuit) or intrinsic (structural de- 
fects inside the material) . Relaxation involves energy ex- 
change between the qubit and its environment and, as 
such, is gencrically sensitive to the low frequency part 
of the environmental spectrum. Relaxation also leads to 
decoherence defined as the decay of off diagonal matrix 
element in the qubit's eigenbasis. 

But dephasing can also occur without energy exchange 
(pure dephasing). Early circuit-QED experiments used 
a Cooper pair box. The coherence properties of these 
devices arc limited by pure dephasing induced by low 
frequency noise. It arises from fluctuating charges in the 
insulating amorphous AI2O3 used to make the Joseph- 
son junctions [23. It is known that the effect of this 
low frequency noise cannot be described within a Marko- 
vian framework (26[. A possible escape to this problem 
is to operate the qubit at special point where sensitivity 
to voltage fluctuations is at second order instead of first 
order [27|. At this working point, treating the effect of 
low frequency noise requires going beyond the Markovian 
approximation [28|, H^] . Fortunately recent circuit QED 
experiments use another type of qubit, called the "trans- 
mon" which has very low sensitivity to voltage fluctua- 
tions (sweet spot everywhere), can be operated as a two 
level system and the low frequency noise [23| . Therefore, 
having in mind future experiments performed with trans- 
mons, qubit relaxation and dephasing will be treated 
within the Markovian approximation in the present pa- 
per. 

Under this hypothesis, the evolution of the qubit + 
cavity system reduced density operator is described by a 
master equation of the form: 

l—iwrf+j^w do) 

where H denotes the Jaynes-Cummings Hamiltonian ([1]) 
and Cj (p) are the Lindbladian superoperators associated 
with each markovian dissipation channel. 

The relaxation Lindbladian operator is given by: 

A-(p)=7i (*-.p.v + - \ {^Y~>p}) ■ ( U ) 

where 7! represents the relaxation rate of the isolated 
qubit. Markovian pure dephasing corresponds to the di- 
fusion of the relative phase between states |±). The pure 
dephasing Lindbladian operator is then given by: 

C r (p) = ^(a*.p.a*-p) . (12) 

where 7^ is the pure dephasing rate of the qubit defined 
as the decaying rate of the (+\p(t)\—) off diagonal matrix 
element for a standalone qubit. 
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2. Resonator relaxation 

In the case of the circuit-QED devices, relaxation 
mainly comes from capacitives losses at the ends of 
the resonator. The temperature dcpcndance can be ex- 
plained by considering that the inverse quality factor Q _1 
is a sum of two contributions. The first contribution rep- 
resents thermal breaking of Cooper-pairs in the super- 
conductor and it scales as exp (T c /T). The other contri- 
bution has a much weaker temperature dependance and 
is probably associated with intrinsic dissipation mecha- 
nisms such as dielectric losses or magnetic vortices [2lj . 
This term is responsible for the low temperature value 
of quality factor. Note that the current limitation in 
the quality factor of microwave cavities used in Rydberg 
atom experiments comes from diffraction losses, surface 
rugosity and residual resistance of the superconducting 
mirrors coming from defects, impurities and magnetic 
vortices [30( ] . 

The corresponding Lindbladian is then: 

C c (p) = k ^a. p. a} — — (at a. p + p. a) a^j . (13) 
where k denotes the total relaxation rate of the oscillator. 



corresponding to a i rotation around the z axis, Lji = 
y^fi<J~ and L c — y/na. 

Between these jumps, the evolution is described by an 
effective Hamiltonian that describes both its intrinsic dy- 
namics and the acquisition of information arising from 
the fact that no quanta has been detected: 

ffe ff =#-fX>]V (15) 

3 

The reduced density matrix is then recovered by averag- 
ing over the set of stochastic trajectories associated with 
a large set of quantum jumps sequences. The weight of a 
given trajectory can be directly related to the dates and 
types of the various quantum jumps. The corresponding 
Linbladian is then given by: 

C j (j>)=L j .p.L\-±{L\L ii p}. (16) 

This method proves to be very convenient numerically 
since the number of variables involved is of the order of 
the dimension d of the system's Hilbert space whereas it 
scales as d 2 in the master equation approach. 



IV. DISSIPATIVE DYNAMICS 



B. Pointer states dynamics 



A. The stochastic wave function method 



In principle, cq. (|10p can be solved numerically in order 
to obtain the quantum dynamics. However, an analyti- 
cal ansatz for the reduced density matrix can be found 
within the mesoscopic approximation using the the quan- 
tum jump approach 12| to the dissipative dynamics of 
the atoms + cavity system. As we shall see extensively in 
this paper, it provides a deep and useful insight into the 
full dissipative dynamics of the cQED system and will 
enable use to find simple analytical results for the Rabi 
oscillation signals. 

The quantum jump method provides a solution to the 
master equation (|10p by assuming that the environment 
of the system is continuously monitored so that any emis- 
sion or absorption of quanta by the system can be as- 
signed a precise, although stochastic date. Each time 
such an event occurs, the system undergoes a quantum 
jump: 



m + )) = 



( L j L 3)mt-)) 



(14) 



The probability rates at a given time t for the vari- 
ous quantum jumps are directly obtained as averages 



where the Lj denote the quantum jump 



~j^3/m-)) 

operator of type j. Here, the quantuum jump operators 
associated to pure dephasing, qubit relaxation and cav- 
ity relaxation are respectively equal to: L v = y // "/ lp /2ia z 



A general stochastic dynamics tends to produce arbi- 
trary mixtures of states. However, it takes a very sim- 
ple form for states that are preserved (up to a phase) 
between and during quantum jumps. These so called 
pointer states entangle minimally with the system's en- 
vironment. In some cases, they provide a basis of the 
system's Hilbert space, thus leading to a simple descrip- 
tion of the stochastic dynamics which we call a pointer 
state dynamics. 

Since coherent states are pointer states with respect to 
cavity losses, such a description is directly relevant for 
the describing the effect of cavity relaxation on cavity 
QED systems in both dispersive and resonant regime. 



1. General results 

We assume that each eigenvalue of the quantum jump 
operator L is non degenerate and choose a fixed section 
A m | A) {L\X) = A | A}) of the corresponding vector bun- 
dle over L's spectrum. Pointer states |A) remain pointer 
states in the evolution between quantum jumps as soon as 
[L,H] = h(L) and [L,L^L] = f(L) (sufficient condition). 
In this case, evolution beween quantum jumps of |A) 



produces a single quantum trajectory t 
where: 



d\(t) 



,i0(\,t) 



dt -£*W*))-i/(A(*)) 



|A(t)> 



(17) 
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with initial condition A(0) = A. The phase 0(X,t) con- 
tains an Hamiltonian and a Berry phase contribution: 



6 t [X] = -lj\x(r)\H\X(r))dr 



(18) 



3((A T |A(r)»dT 



(19) 



Then, starting from a single pointer state |A), the 
stochastic dynamics taking into account all possible se- 
quences of quantum jumps occuring at times < t\ < 
. . . < t p < t produces a single trajectory 1 1— > A(i) but en- 
codes the sequence of quantum jumps in a (t\ , . . . , t p ) de- 
pendent phase. The resulting state at time t is of the form 
e t9(ti, ...,*„) \X{t)) and therefore, any initial state which is 
linear combination of pointer states will experience deco- 
herence because of the averaging of these random phases. 

Summing over all quantum jumps sequences leads to 
the system's reduced density operator obtained from an 
initial state c\ |A): 

P (t)= J2 c A+ cL^[A + ,A_] e ^[ A +l- e '[ A -])|A + (t))(A 
(A+.A-) 

(20) 

where the decoherence functional is given by: 

^[A + A_] = e -i/o |A + -A_| 2 (r)dr e47 / f Q(A + Al)(r))dr 

(21) 

As expected, this expression coincides with the accumu- 
lated decoherence of a pair of coherent states driven along 
trajectories 1 i— > X±(t). 



substituting a±(r) into (|21|) . As we shall see in more de- 
tails in section HVD 2i this is exactly the dynamics solved 
in Ref. [l3| using a different method. 



In the case of a qubit resonantly coupled to the cav- 
ity, generalized Gea-Banacloche states are approxi- 
mate pointer states with respect to cavity losses. Let 
us first consider an initial state of the form \ip q b) <& |a) 
where \ip q b) is the qubit initial state and where the 
parameter a of the initial coherent field is taken real. 
Cavity losses will then produce a qubit + cavity re- 
duced density operator involving two generalized Gea- 
Banacloche states \ty±(t)). The decoherence coefficient 
in front of the (t)\ is obtained by substitut- 

ing a±(t) = a e™*/ 4 ^ in (f2"Tj) . Introducing one or 
several 7r-pulses (single echo or bang-bang control) only 
changes the trajectories of the Gea-Banacloche param- 
eters a±(t) in the Fresnel plane but eq. (|21[) remains 
valid. To deal with more general initial condition such as 
■0qb) <& (|«i) + \ a 2)), one just has to decompose the ini- 
- Wtial state as a sum of rotated generalized Gea-Banacloche 
states where the X direction is replaced by the directions 
given by the phases of the coherent states present in the 
initial condition. 



So far, only the effect of cavity losses has been dis- 
cussed. How do other sources of dissipation affect this 
picture in the dispersive and the resonant regimes ? 



2. Application to cQED systems 

In the dispersive regime of cQED, the effective dis- 
sipative dynamics is described by the effective Hamil- 
tonian ([2]). In the strong coupling regime g <C k, the 
effective quantum jump operator associated with cavity 
losses can still be approximated by L = ^/Ha. There- 
fore, in this scheme, the pointer states with respect to 
cavity losses are states of the form |±) ® \a). Since any 
state can be expanded as a linear combination of these 
pointer states, the dispersive regime of cQED realizes the 
above describes pointer state dynamics. Their respective 
coherent state parameters evolve according to: 

^ = -iw±a±(t)-|a±(t)-ie(t) (22) 

where lo± = loq ± x ' IS the ac-Stark shifted cavity fre- 
quency and e(t) is a driving of the cavity. Starting from 
an initial state of the form IV'qb) <S> \a) leads to a qubit + 
cavity reduced density operator of the form PO)) involv- 
ing two pointer states |+) <X> \a + (t)) and |— ) <g> |a_(t)) 
where a±(r) denote solutions of (|22p with initial condi- 
tion a±(0) = a. The decoherence functional induced by 
cavity losses and relative to these states is obtained by 



To answer this question, let us first remark that, in 
the dispersive regime, pointer states with respect to cav- 
ity losses are also pointer states with respect to pure de- 
phasing of the qubit. Qubit relaxation sends a pointer 
state of the form |+) (g) \a) on a pointer state |— ) <g> |a). 
Therefore, in the dispersive regime, a pointer state is sent 
on another pointer state. But the situation is more in- 
volved in the resonant or quasi-resonant regimes. How- 
ever, we shall see that quasi pointer states are sent on 
linear combinations of quasi pointer states. Thus, the 
complete stochastic dynamics can still be described in 
terms of generalized Gea-Banacloche states © and for 
this reason, simple analytical results can still be obtained. 

The simpler case of pure dephasing will be considered 
first and we will then turn to the more complicated case of 
qubit relaxation. Explicit analytical results for the spon- 
taneous and induced Rabi oscillation signals are given 
in section IIVEI For the sake of simplicity, the stochas- 
tic dynamics will be discussed using the factorized form 
d3J of generalized Gea-Banacloche states, assuming that 
their coherence is not broken by the dissipative processes 
under consideration. This assertion will be justified by 
a more precise discussion of the dynamics postponed in 
appendix |Bj 



C. Pure dephasing 

1. Resonant regime 

First of all, let us note that the corresponding quantum 
jump operator has no effect between quantum jumps. 
Moreover, the rate of pure dephasing jumps is constant 
in time, independent of the state and equal to 7^/2. 
Next, each quantum jump acts as a 1 rotation around 
the z axis, exactly like an echo pulse. As such, it sends 
the atomic polarization \D±(8)) defined in cq. ([5]) on 
\D T (—6)) not affecting the field component ([6|). Af- 
ter p jumps occuring at times < t\ < ... < t p < 
t, an initial state \D m (0)) <B> \tp m (0)) has turned into 
e- imeVW \D m ,(9 t )) <E> \ip m >(0 t )) where m' = (-1)% and 
8t = 8{ti, . . . , t p ; t) is defined as 

6(t 1 ,...,t p ;t)=gf2(-l) j (tj+i-tj) (23) 

3=0 

where to = and t p+ \ = t. 

A typical trajectory T of the associated phase in Fres- 
ncl plane is depicted on figure [IJ with slope changing of 
sign at each quantum jump. The relative phase between 
the two quasi-coherent field components follows a random 
trajectory directly related to the integral of a telegraphic 
noise: 

(A0)(t) = -4= f X(r)dr (24) 
2y/n Jo 

where X(t) G { — 1, 1} starts with X(0) = 1 and jumps 
each time there is a quantum jump. The probability dis- 
tribution of waiting times r between quantum jumps is 
exponential i^(t) = ^ e _7vT / 2 . The characteristic func- 
tion of the probability distribution for ( A</>) (i) is nothing 
but the decoherence coefficient of a qubit in a longitudi- 
nal telegraphic noise [3l[. Its behavior depends on the 
dimensionless parameter r\ = g/-fip \ffi which is equal to 
2ix / N v where N v = 2ir\/E'y ip /g is the average number 
of quantum jump events occuring before the date of the 
first expected Rabi oscillation revival. 

For 77 3> 1, only a single event is necessary to spread the 
phase (A</>)(i) over 2ir. In this regime, the evolution of 
the cQED system up to the first revival is dominated by 
quantum trajectories having no pure dephasing quantum 
jumps. Their total weight is e~ lift l 2 ~ 1. Trajectories 
with one quantum jump or more have a total weight 1 — 
g-7.pt/ 2 <g; 1 anc j can therefore be treated perturbatively. 
For fixed g and -f v , this " weak dephasing limit" is realized 
for n <C {g /lip) 2 ■ In practice, this is the regime of interest 
for present and forthcoming experiments since g/j v is 
above 40. 

For i)< 1, a large number of events are necessary to 
spread significantly the random relative phase (A</>)(i). 
In this case, rapid dephasing of the qubit prevents the ob- 
servation of an atom + cavity entangled state. As shown 



7 




tx t 2 



FIG. 1: Typical trajectory of the phase cj> of the quasi-coherent 
component state in the Fresnel plane showing two quantum 
jumps at time t\ and ta • The effect of a pure dephasing quan- 
tum jump is to reverse the direction of motion of the quasi- 
coherent component in the Fresnel plane. Inset recalls the 
definition of </>. 



in appendix [B] the incoherent qubit ends up breaking 
the coherence of the coherent field, selecting Fock states 
as pointer states. This is the "strong dephasing limit" 
which, as already stressed, can only be reached for very 
large photon numbers in the strong coupling regime of 
cQED. 

Nevertheless, in the strong dephasing limit, even if 
the spreading of the field phase takes place over a time 
scale 2j v n/g 2 ^> the Rabi oscillation signal will 
be strongly damped in a much shorter time. Indeed 
for t < 2ir\/E/g, although the Fresnel angle of Gea- 
Banacloche states are weakly dispersed around zero and 
thus have a rather strong overlap for a typical given 
quantum trajectory, the classical Rabi oscillation phase 
e 2in<t> t avera g CS to zero over a time scale of the order of 
27" . This means that the qubit dephasing time is any- 
way the upper limit to the Rabi oscillation visibility in 
the mesoscopic regime. Nevertheless, for qubits in high 
quality resonators but rather poor dephasing properties 
g/'Ytp < 5, the crossover from weak to strong dephasing 
regimes might be observed by performing a tomography 
of the field state. 

Finally, the result of this analysis is that, in the strong 
coupling regime of cQED, the main contribution to the 
Rabi oscillation signal comes from quantum histories 
without any pure dephasing quantum jump. As far as 
we are interested by this signal, the effect of qubit de- 
phasing can be accounted for through a supplementary 
decoherence factor on top of the decoherence factor as- 
sociated with cavity losses: 

^l(t)~e-T^ 2 . (25) 



2. Dispersive regime 

In the dispersive regime, the analysis is much sim- 
pler since pointer states with respect to cavity losses are 
also pointer states with respect to pure dephasing of the 
qubit. Therefore, starting with an initial state of the form 
iV'qb) <8> | a), the only effect of pure dephasing is to mul- 
tiply the decoherence functional associated with cavity 
losses by: 

jf}(t)=e-^ t . (26) 

D. Qubit relaxation 

1. Resonant regime 

For a standalone qubit, the non hermitian part of the 
Hamiltonian tends to bring the state of the qubit towards 
the |— ) state, thus reflecting the acquisition of informa- 
tion associated with the absence of relaxation quantum 
jump. But at resonance, in the presence of a meso- 
scopic state, the pumping of the qubit by the cavity 
photons alters the dynamics between quantum jumps. 
As discussed in appendix \^ in the 71 C 3 regime, 
the dynamics is dominated by the strong atoms + cav- 
ity coupling. Contrarily to the case of a the dispersive 
regime where no photon can excite the qubit after re- 
laxation, several relaxation jump can occur. The statis- 
tics of relaxation quantum jumps is indeed described by 
a renewal process with effective waiting time distribu- 
tion Vi(r) = (7i/2)e" 7lT / 2 . The effective rate 71/2 fol- 
lows from the localization of generalized Gea-Banacloche 
states in the equatorial plane of the Bloch sphere. 

When the first relaxation quantum jump occurs at 
time ti, the qubit gets projected on the |— ) state leav- 
ing the cavity in a superposition of quasi-coherent states 
IV'i^i))- The subsequent Rabi oscillation signal then 
shows rapid oscillations of frequency g\ffi which corre- 
spond to the superposition of the Rabi oscillation signals 
associated with the evolution of |— ) ® \ip±(ti))- There- 
fore, the first relaxation jump leaves us with four quasi- 
coherent components in the Fresnel plane as shown in fig. 

HI 

o |**(0)) = e —^- - |**(-0)» (27) 

The Rabi oscillations starting right after the first jump 
collapse in time g~ x due to the splitting of the pairs of 
counter-rotating quasi components in the Fresnel plane. 
Revivals then occur when two of them recombinc. 

Assuming that gt\j2\ffi < tt, the first recombination 
occurs at time 2ti and involves components (b) and (c) on 
fig. [2]). Components (a) and (b) will recombine at time 
t\ + Airy/E/g as well as components (c) and (d). But in 
these three cases, the rapid oscillations are simply washed 
out by the averaging over t\ as in the pure dephasing 
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FIG. 2: Quasi-coherent components generated after one relax- 
ation jump pictured in the Fresnel plane. Field quasi-coherent 
states are represented as uncertainty disks at the tip of their 
classical amplitudes. Components components (b) and (d) 
move clockwise whereas components (a) and (c) move coun- 
terclockwise. 



case. Only the recombination of components (a) and (d) 
will not be averaged to zero since it takes place at time 
Air^/K/g, independent of t\. 

A careful discussion of of all contributions is presented 
in appendix [Bl It shows that keeping only the first term 
in the r.h.s. of (|27|) corresponds to retaining only tra- 
jectories for rapid oscillations survive the averaging over 
dates of relaxation jumps. Therefore, we are brought 
back to the case of a single pair of trajectories within 
the Fresnel plane. But with each relaxation jump comes 
not only a pure phase as in eq. (38) of [10( but a co- 
efficient e -1 ^'™. 8 / '\/2. Following [loj |. these coefficients 
can be resummed thus leading to the decoherence coeffi- 
cient associated with relaxation in front of the coherence 
|*?(f)X*5(t)|: 

4 fl) (*)=exp(f |(i e "*M-l)dr) (28) 

where (A</>)(r) = gr/2^/E denotes the angular separation 
in the Fresnel plane of the quasi-component components 
of|*£(r)>. 

As a final comment of this discussion, let us point out 
that this decoherence coefficient alone is not sufficient 
to keep track of the complete state of the atoms + cav- 
ity system. More precise results could be obtained by 
perfoming a pertubative expansion in 71 / g keeping track 
of all contributions arising from (|27|) . Cavity losses and 
dephasing jumps can then be taken into account. The 
resulting formula are appropriate in the weak relaxation 
regime but rapidely require numerical evaluation. 
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2. Dispersive regime 

Being off resonance, the cavity cannot send back the 
qubit to the |+) state after a relaxation quantum jump. 
As stressed before, since the states |±) (g) \a) are pointer 
states with respect to cavity relaxation, the relaxation 
dynamics is most conveniently studied in term of these 
states. 

The state |— ) <g> \a) are obviously left invariant by the 
relaxation dynamics since a~ . \— ) <£> \a) = 0. The state 
|+)®|a) is sent on |— )<g>|a) by the relaxation jump. The 
probability for such a jump to occur between t and t + dt 
is 71 dt, independent of a. With this in mind, an exact 
solution for the dynamics in the dispersive regime, taking 
into account all dissipative processes (relaxation of the 
cavity, qubit relaxation and dephasing) can now be given 
for an initial state of the form (^4 + |+) + }) <S> 

Following Ref . [HI , the qubit + cavity reduced density 
operator can be decomposed with respect to the qubit 
state p(t) = Y,ee' Pe,e'(*) ® | e) (e' | . Note that p_ + = 

p\ and Tr((p ++ + p ){t)) = 1. First of all, diagonal 

terms p e ,e{t) are not affected by pure dephasing. Next, 
the p+ t +(t) operator only contains contributions coming 
from trajectories without any relaxation quantum jump 
whereas p_._(t) is fed from trajectories that have had no 
quantum jump and trajectories having a single jump at 
time < r < t. This leads to: 



E. Rabi oscillation signals 

1. General form of the signal 

Let us now use these results by computing the Rabi 
oscillation signal and its enveloppe assuming that the 
average photon number in the cavity remains equal to 
n. The result for the Rabi oscillation signal P(t) defined 
as the probability of finding the qubit in the |+) state is 
given by: 



Pit) = \ 



1 + 8? 



(32) 

where TZ^ (t) = (ip~(t)\ip + (t)) measures the overlap be- 
tween the two Gea-Banacloche quasi coherent component 

of the field. Dccohcrcncc is contained in (t). Given 

these coefficients, the upper and lower enveloppcs P± of 
the Rabi oscillation signals are given, in the case of a 
single qubit initially in the excited state, by fTo|: 

P±(t) = ±(l±\K + -(t)F+-(t)\) . (33) 

Explicit expressions for the overlap and decoherence co- 
efficients will now be given for the case of free evolution 
of the cQED system (spontaneous revivals) and for the 
case of an echo experiment (induced revivals). 



p++(t) = \A + \ 2 e-^\a + (t))(a + (t)\ 
p__(t) = \A.\ 2 \a.(t))(a.(t)\ 

-7l* 



(29) 



7i e 



(t,T))(a-(t,T)\dr (30) 



2. Spontaneous revivals 

In the case of a free evolution (spontaneous revivals), 
the overlap factor (i) is given by: 



where a±(t) denotes the solution to ([22]) with corre- 
sponding sign and initial condition a± (0) = a and 
a.- (t, t) denotes the solution of l|22p with + sign from 
to t and then solution to (f2"2"| with — sign. 

Because coherence operators p^ (t) and p |_ (t) in- 
volves a |— ) state, they only receive contributions from 
trajectories without any quantum jump. Taking into ac- 
count pure dephasing leads to: 



x e 



i(<p + -(p_)(t) 



«+(*)><«-(*)! 



(31) 



These formula arc nothing but generalizations of eqs. 
(5.14) to (5.18) of [HI that take into account relaxation 
of the qubits and for a general driving of the cavity. 

Our derivation, not relying on the P-function formal- 
ism, provides a simple view of the underlying quantum 
dynamics showing the importance of pointer states with 
respect to to cavity losses and qubit pure dephasing. Fol- 
lowing the same line of reasoning, the reduced density op- 
erator for several qubits coupled to a driven cavity all in 
the dispersive regime can be computed although it leads 
to more complicated expressions. 



Tl + -(t) = e~ n e 



n igty/n 



+ °° 



E 

fe=0 



-igtvE+1 



k\ 



(34) 



Taking into account all decoherence sources, the decoher- 
ence factor JF_| (t) = e - (i ( t )+ 40 ( t ) i s given by: 



d(t) 



nn H — — I t 



(«n+^) sin(^t) (35) 



where <p t — gtj1\fn and 



Q(t) = -—- sm 



(36) 



Decoherence is mainly dominated by an exponential de- 
cay with rate T = Kn+("f lp +"/i)/2. The appearance of nn 
reflects the photon emission rate enhancement by stimu- 
lated emission. The 1/2 reduction of the relaxation rate 
comes from the averaging over classical Rabi oscillations 
of the qubit relaxation rate. To explain the 1/2 reduc- 
tion of the dephasing let us note that, for a standalone 
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qubit, contributions from all quantum history pile up to 
give (+|p(t)|+) = (+|j°(0)|+) e -7 ^'. On the contrary, for 
a qubit coupled to a mcsoscopic field in a cavity, only 
quantum histories without any quantum jump before the 
measurement time contribute to the Rabi oscillation sig- 
nal. 



3. Induced revivals 

In the echo experiment with echo pulse performed at 
time tir, the Rabi oscillation signal at time t < is 
determined by eqs. d34j), (35|) and (36|). For t < in, 



the overlap factor is given by (t w ,t) = (2t n — t) 
where the r.h.s. involves the free evolution overlap factor 
since the echo pulse reverses the dynamics of the qubit 
+ cavity state |9j. 

The decoherence factor = e -d{U,t)+ie(t«,t) 

for t > t n is now given in terms of: 



(ku + { sm - <f>t) - 2 sin (<f> w )} 



(7». + 7i)) t 



(37) 



where <f> n = gt 7r /2\/n. Its phase is the sum of a contri- 
bution due to cavity relaxation (eq. (51) of [lj|) and a 
contribution due to qubit relaxation: 



7i 1 



2 sin 2 



(38) 

Note that the dephasing exponential factor e t^'/ 2 is the 
same than in the free evolution since the echo pulse is not 
able to reverse the effect of a Markovian qubit dephasing. 

Besides the overlap coefficient TZ-\ (t n , t), the echo pulse 

mainly affects the contribution arising from the accumu- 
lation of slow phases either associated to cavity relaxation 
or to qubit relaxation. 



V. DISCUSSION OF THE RESULTS 

A. Methods and parameters 

The Yale group reports of values of g/2n between 5.8 
and 100 MHz [32||. Low relaxation and dephasing rates 
have recently been measured [ij at the magic point n g = 
1/2: k/2tt ~ 0.6 MHz, j 1 /2n ~ 0.02 MHz and 7 v /2tt ~ 
0.31 Mhz but for a low value of the coupling [38| : g/2ir ~ 
11.6 MHz (see row Circuit QED (1)). In a more recent 
experiment 0], a coupling of g/2n ~ 210 MHz has been 
obtained with dissipation parameters k/2it ~ 0.25 MHz, 
7i/27r ~ 1.8 MHz and ^ v /2ir ~ 1 Mhz (see row Circuit 
QED (2) in table H}. The last row of table Q] present 
choices of dissipation parameters compatible with near 
future improvements of circuit QED devices. 



Type of device 




5/7i 


g/iv 


Rydberg atoms (1) 


310 


10230 


00 


Rydberg atoms (2) 


4300 


10230 


00 


Circuit QED (1) 


19.4 


580 


40 


Circuit QED (2) 


840 


106 


215 


Circuit QED (3) 


1400 


2000 


2000 



TABLE I: Physical parameters for cQED experiments per- 
formed with Rydberg atoms (LKB group) or with supercon- 
ducting circuits (Yale group) . Experimentallyrealized config- 
urations appear on rows Rydberg atoms (1) [8J|, Circuit QED 
(2) [| and Circuit QED (2) Q. 



Analytical results presented in section IIVEI have been 
compared to quantum Monte-Carlo simulations of the 
qubit + cavity evolution. For these simulations, the 
Adams-Blashford scheme of order four has been used to 
compute the evolution under the non-hermitian Hamil- 
tonian between quantum jumps. All simulation results 
presented here have been obtained using an average over 
10000 trajectories. 

Results for the free evolution of the qubit + cavity 
system will be first presented in section IV Bl and for the 
echo experiment in scction fV CI Consequences from these 
results are drawn in section |V Dl 



B. Free evolution 

As a first step, it is interesting to see how our ana- 
lytical model provides accurate predictions in the pres- 
ence of qubit relaxation and pure dephasing on top of 
cavity relaxation. Figure [3] presents a comparison of 
the analytical enveloppes for the Rabi oscillation signal 
S z (t) = P m =in{t) — 1/2 in the presence of (a) cavity 
relaxation alone, (b) cavity and qubit relaxation but no 
pure dephasing, (c) cavity relaxation and pure dephasing 
and finally (d) including all dissipative mechanisms. The 
dissipative parameters used for this comparison corre- 
spond to the Circuit QED (2) row of tablcQ] g/n = 840, 
17/71 = 106 and g/j v = 210. With these parameters, the 
analytical model correctly describes the enveloppe of the 
Rabi oscillation signal. 

Figure 0] compares the Rabi oscillation signals for the 
experimental parameters corresponding to Circuit QED 
lines in table HI In cases (a) to (c), our anaytical model 
predicts the upper and lower enveloppes for the Rabi os- 
cillation signal with good precision. But in case (d), only 
the initial collapse is correctly accounted for by the ana- 
lytical enveloppes. Our Monte-Carlo simulation shows a 
relaxation of S z (t) towards —1/2 at longer times which is 
not described by our analytical model. This result is not 
surprising since for g/n = 19.4, at t = Mr, the average 
photon number in the cavity should have decayed from 
10 to 2 because of cavity relaxation alone. Therefore, at 
t ~ 5tR, the qubit + cavity system is already out of the 
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FIG. 3: (Color online) Influence of the various dissipation 
mechanisms on the spontaneous revivals for a single qubit in 
the presence of a n = 10 photon field. The graph on the 
left depicts the analytical enveloppes for the Rabi oscillation 
signal S z (t) — P m= i/2(i) — 1/2 as a function of t/tn for no 
dissipation (black lines), (a) only cavity dissipation g/n — 
840, (b) cavity and qubit relaxation g/"fi = 106, (c) cavity 
relaxation and qubit pure dephasing g/^ v = 210 and (d) all 
dissipation mechanisms present. The right part of the figure 
presents the associated Rabi oscillation signals obtained from 
quantum Monte-Carlo simulations (plain line) as well as the 
corresponding analytical enveloppes in cases (a) to (d). 



FIG. 4: (Color online) Spontaneous revivals for a single qubit 
in the presence of a n = 10 photon field for various dissi- 
pative parameters. The graph on the left depicts the an- 
alytical enveloppes for the Rabi oscillation signal S z (t) = 
P m -i/2(t) — 1/2 as a function oit/tn for (a) no dissipation, (b) 
Circuit QED (3) parameters, (c) Circuit QED (2) parameters 
and (d) Circuit QED (1) parameters (see table UJ) . The right 
part of the figure presents the associated Rabi oscillation sig- 
nals obtained from quantum Monte-Carlo simulations (black 
lines) as well as the corresponding analytical enveloppes in 
cases (a) to (d). 



mesoscopic regime and our model is not expected to 1 
valid in this regime. 

Finally, these numerical and analytical results sho 
that observing spontaneous Rabi oscillation revivals : 
circuit QED is not possible with the Circuit QED (: 
parameters of table |TJ Only the recent improvements c 
the value of g/n (Circuit QED (2) parameters) open tl 
possibility for observing this phenomenon. 

The influence of photon number can be visualize 
through the effective decoherence coefficient C(t, n) 

(t)\ as a function of t/tfi = gt/2ir and n. It co] 

tains the effect of the environment on the contrast 
Rabi oscillation revivals. Contour plots C(i, n) = 0. 
0.8, 0.7, 0.5 and 0.2 are shown on figure[5]for the Circu 
QED parameter sets of table HI The fourth graph sho'v 
these coutours for the case of a Rydberg atom expcrimci 
assuming a cavity damping time of 14 ms. 

At high photon number, the various coutour curves 
tend to become vertical. This is easily understood by 
noticing that, at high photon number, the cavity deco- 
herence coefficient becomes independant of n and domi- 
nated the decoherence process. In this regime, C K (t) ~ 
cxp (— 2-7T 3 (g/n) (t/tn) 3 ) and the decoherence time scales 
typically as irtu{2g / n) 1 ^ . Figure [5] shows that a co- 
herent mesoscopic dynamics could be observed in really 
good conditions in cases (c) and (d): C(t,n) > 0.8 at 
t = t.R\ffi corresponding to the generation of a meso- 
scopic field Schrodingcr cat state involving up to 15 pho- 
tons. 
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FIG. 5: (Color online) Coutour plots of the effective decoher- 
ence coefficient (t)| as a function of t/tn and n (n > 5 to 

ensure validity of the mesoscopic approximation). The first 
three graphs represent the case of circuit QED devices: (a) 
Circuit QED (1) parameters, Circuit QED (2) parameters and 
(c) Circuit QED (3) parameters. Graph (d) depicts the cor- 
responding information for Rydberg atom experiments with 
g/ti = 4310 corresponding to a cavity with 14 ms damping 
time. In this case we extended the n range up to 50 pho- 
tons. On each graph, the dashed line corresponds to the time 
of Schrodinger cat preparation t = tR^ffi and the dot-dashed 
line correspond to the maximum of the first spontaneous Rabi 
oscillation revival t = 2tR\ffi. 
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FIG. 6: (Color online) Induced revivals signal S z (t) for oi 
qubit in the presence of a n — 10 photon field for varioi 
dissipative paramerters as a function of t/tn. Graphs (a) 
(c) depict Rabi oscillation signals obtained from a quantum 
Monte-Carlo simulation for t^ = 3ta (back curve). Analyt- 
ical enveloppes showing the contrast at the induced revival 

(t/2, t)\ taking into account all sources of dissipation 

(full blue lines) or retaining only cavity relaxation (dashed 
blue lines) are also displayed. Graph (d) depicts the contrast 
of the induced revival as a function of t/tn. — 2t n /tR using 
parameters of table [J The Circuit QED (n) parameter set 
corresponds to curve Cn (n = 1 to 3) and the Rydberg atom 
(m) set parameter corresponds to curve Rm (m = 1,2). 




FIG. 7: (Color online) Coutour plots of the effective echo de- 
coherence coefficient {t/2, t)\ as a function of t/tn and 
n (n > 5 to ensure validity of the mesoscopic approxima- 
tion). The first three graphs represent the case of circuit 
QED devices: (a) Circuit QED (1) parameters, Circuit QED 
(2) parameters and (c) Circuit QED (3) parameters. Graph 
(d) depicts the corresponding information for Rydberg atom 
experiments with g/n = 4310 corresponding to a cavity with 
14 ms damping time. In this case we extended the n range 
up to 50 photons. 



C. Echo experiment 

First of all, we have considered a circuit QED echo ex- 
periment obtained with 10 photons initially and dissipa- 
tive parameters g/n = 840, 3/71 = 106 and g/j v = 210 
(Circuit QED (2) of table H}. Figure [6] shows the echo 
signals for t/tn = 3 as well as the analytical enveloppes 
obtained for t = 2t„ when the two quasi-coherent com- 
ponents of the field recombine. 

Because the echo experiment reverses exactly the 
dyamics of the qubit + cavity system, the reduction 
of the induced Rabi oscillation revivals contrast is ex- 
actly given by the environmental decoherence: C e (t, n) = 
\J r + -(t/2,t)\. Figure [7| shows coutour plots for C e (t,n) 
as a function ot t/tn and n using the sets of parame- 
ters given in table [T] for circuit QED devices. As in the 
previous paragraph, the corresponding plot for Rydberg 
atom experiments with cavity damping time of 14 ms are 
given. 



D. Discussion of the results 

Within the context of circuit QED experiments, our re- 
sults suggest that tests of quantum coherence in the res- 
onant mesoscopic regime could be performed using the 
latest generation of circuit QED devices. Echo exper- 
iments could be performed with up to 10 photons and 
echo times up to five vacuum Rabi periods. A more pre- 



cise insight into the field dynamics could be obtained 
by performing a field tomography giving access to either 
the Q or the Wigner function of the field. The Q func- 
tion probes the phase distribution of the field and its 
experimental determination would test the splitting of 
the field into two quasi-coherent components with oppo- 
site phases. Observing the induced revival then provides 
an experimental proof of the generation of a mesoscopic 
superposition superposition in the cavity Q . 

On the other hand, Rydberg atom devices [33| are 
characterized by very long atomic relaxation time (33 ms) 
and the absence of atomic dephasing. The main limi- 
tation comes from cavity relaxation and time of flight. 
Up to 2005, the LKB team was able to work with a 
~ 1 ms cavity relaxation time corresponding to 
g/n ~ 310 [|[ (see row Rydberg atoms (1) in table [J). 
Cavity damping times up to 14 ms have been consid- 
ered in our recent publication [Io| an are recalled in row 
Rydberg atoms (2) of table |TJ Recently, the LKB team 
reported on the realization of a ultra-high finesse cavity 
reaching k _1 ~ 10 Hz [3(J corresponding to g/n ~ 31000! 
With such performances, the only limitation is the atomic 
time of flight. This is an important limitation for the 
study of spontaneous revivals. In practice, only partial 
spontaneous revivals arising from small atomic ensembles 
could be observed [l(|. It is worth mentioning that cir- 
cuit QED experiments do not suffer from this limitation 
although their dissipative properties arc still behind the 
Rydberg atom ones. 
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VI. CONCLUSION 

Motivated by rapid experimental progresses of circuit 
QED devices, we have studied the dissipative dynamics 
of one qubit resonantly coupled to a mesoscopic field in 
a cavity taking into account cavity relaxation as well as 
qubit relaxation and pure dephasing. The resonant cou- 
pling between the qubit and the cavity produces an en- 
tangled state involving quasi-pointer states with respect 
to cavity relaxation. This very precise fact ultimately jus- 
tifies in full generality the simple analytical model previ- 
ously used to describe the dynamics of Ry dberg atom 
ensembles coupled to a high-Q cavity [1Q|. Although 
quasi-pointer states naturally produced by the resonant 
interaction are not preserved by qubit relaxation and de- 
phasing, we have obtained a simple description of the 
dynamics of the cavity + qubit system in the presence 
of qubit dissipative processes. Simple analytical expres- 
sions for the spontaneous and induced Rabi oscillation 
signals have been obtained and shown and to be in good 
agreement with quantum Monte-Carlo simulations. 

Our study shows that these effects and the genera- 
tion of mesoscopic Schrodinger cat states could be within 
experimental reach in forthcoming circuit-QED experi- 
ments. However, our prediction may still be over opti- 
mistic since we did not attempt at modeling precisely the 
qubit readout. Thus, as done in Ref. [H for the Ramsey 
fringes measurement, the measurement beam can only be 
switched on after the qubit + cavity resonant interaction 
period, once the qubit has been sent off resonance for 
dispersive readout. Obtaining more precise quantitative 
results suitable for describing an actual experiment will 
probably require a more complete modeling of the experi- 
mental run, from preparation to measurement. The same 
remark is also relevant for obtaining precise predictions 
for the result of a cavity state tomography following for 
example a proposed protocol recently proposed in Ref. 
[3^ . Ultimately, we may probably have to rely on quan- 
tum Monte-Carlo simulations. 

Other superconducting circuit for on-chip cavity QED 
have been proposed. In these proposals, the microstrip 
resonator is replaced by another oscillator such as an LC 
circuit or a biased large Josepshon junction in the ap- 
propriate parameter regime. Buisson and Hekking [l4| 
have proposed a circuit that couples a current biased dc- 
SQUID to a Cooper pair box. In principle, values of 
17/277 of the order of 100 to 200 MHz are within reach 
with such a device. However, modeling its dissipative 
dynamics requires taking into account the effect of low 
frequency noises which act on the qubit through low fre- 
quency dephasing and on the resonator through slow fluc- 
tuations of its resonant frequency. These effects are easy 
to account for using our formalism and a quasi-static ap- 
proximation. But the resonator is also coupled to high 
frequency current noise through the SQUID [35|. Work 
is in progress to describe the effect of this coupling on 
the Buisson-Hckking cavity-QED circuit. 
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APPENDIX A: NON UNITARY EVOLUTION 
BETWEEN QUANTUM JUMPS 

The non hcrmitian Hamiltonian describing the evolu- 
tion between two quantum jumps is: 

rtr = /r^ c -f I^-fiv-^ (Ai) 

Obviously, the pure dephasing part leads to an trivial 
overall e -7 ^'/ 2 factor. In the mesoscopic regime, the re- 
sulting Hamiltonian can then be rewritten in terms of the 
su(2) generators J a suitable for describin g th e dynamics 
of generalized Gea-Banacloche states (|3]) (lOj : 

TT'tfeff = g^Jp+lJ x - 

-'-Z^ J' -{<* + «)■ (A2) 

This Hamiltonian is of the same form as the one discussed 
in appendix D of [id] . The resulting dynamics can then 
be inferred straightforwardly. In [id ], since g 3> k, the 
backaction of cavity relaxation on the qubit has been 
neglected and since g <C 71, the term proportional to 
7i l / z shall also be dropped out. Finally, starting from a 
state of the form 

00 -p/2 

|*(o)) = e ^/ 2 ^^ 7T |^(o)) 

and assuming that -C 1, so that the average photon 
number in the cavity remains constant between and t, 
we get: 

00 — v II 

l*nj.(*)> = e-^e-tfe-E n '^^ |* p (0)> . 

p=0 vP- 

(A3) 

The exponentially decaying prefactors are due to the non- 
hermiticity of the Hamiltonian. They reflect the decaying 
weight in time of quantum trajectories without any quan- 
tum jump between and t. The apparently surprising 
appearance of 71 /4 instead of 71 /2 for relaxation quan- 
tum jumps corresponds to time averaging to 1/2 of the 
qubit populations induced by the strong coupling with 
the cavity. Physically, we assume that the number of 
photons is so high and relaxation so weak that the cavity 
mode keeps driving the qubit with the same intensity. 
Note that the resulting effective dynamics correspond to 
an effective Poissonian statistics of relaxation quantum 
jumps. But when the qubit state is close to |— ), the 
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probability rate for such a jump is much smaller than 
when it is close to |+). Indeed deviations from the ef- 
fective Poissonian statistics average out since dissipativc 
processes take place over a much longer time scale than 
rapid oscillations induced by the resonant qubit/cavity 
coupling. 



APPENDIX B: EFFECTS OF DISSIPATION ON 
THE RESONANT CQED SYSTEM 

In this section, the evolution of generalized Gea- 
Banacloche states ^ under the stochastic dynamics in- 
volving qubit relaxation and dephasing is discussed. We 
first derive the effective dynamics resulting from relax- 
ation and pure dephasing quantum jumps. The case of 
pure dephasing can then be treated exactly using renewal 
theory. A formal solution including relaxation in the low 
dissipation, large photon number limit is provided and 
used to extract the dominant contributions to the Rabi 
oscillation signal. This appendix provides a firm basis for 
the disscusions of sections IIVCI and IIVDI 



1. Stochastic dynamics 

a. Evolution of the coherences between generalized 
Gea-Banacloche states 

In order to compute the evolution of the qubit + cavity 
state, it is convenient to look at the evolution of a coher- 
ence between a given pair of generalized Gea-Banacloche 
states with respective indices (m+ , m_ ) . Such a coher- 
ence is initially of the form 



n, 



-(0) 



" E 

P+.P- 



;(p + +p_)/2 



Vp+ 1 p- 



(p-)l 



(Bl) 

As shown in appendix \K\ in the weak dissipation limit 
and in the mesoscopic regime, the evolution between 
quantum jumps is exactly the same as in the dissipa- 
tionless case apart from the normalization prefactor that 
reflects the weight of quantum trajectories not present- 
ing any quantum jump in the considered time interval. 
Moreoever, in the following, we shall restrict ourselved 
to short times (jt <C 1) so that the average photon num- 
ber in the cavity remains constant. The effect of various 
quantum jumps is described by: 



a 



\X%>) = y/p + l/2-m\xM) 
\X { p ] ) = {\X^- 1] ) -\X (p - 1) ))/2 



\x(p) 



\x 



(p) 



(B2) 
(B3) 
(B4) 



In principle, we should follow keep track of all phases 
in each subspace generated by \X± ^). It was proved in 
Ref. [Tol | that for cavity relaxation, this exact dynamics 
could be further simplified so that the resulting stochastic 



dynamics be expressed only in terms of the generalized 
Gea-Banacloche states A similar route will be fol- 
lowed here and recasting eqs. (|B2[) to (|B4|) in terms of 
generalized Gea-Banacloche states using the mesoscopic 
approximation leads to: 



a.\**{9)) 
*-.|*m(0)> 



(B5) 



a i<t>m,e 



V2 

\* x m (-o)) 



(|*?(0)>-|*-HO>)(B6) 



(B7) 



b. Trajectories of generalized Gea-Banacloche states 

A sequence of quantum jumps of the above types will 
generate a tree structure whose vertices are associated 
with the branching expressed in eq. (|B6[) . Vertices asso- 
ciated with photons are located on straight edges of the 
tree and vertices associated with pure dephasing jumps 
are associated with direction reversals (see fig. [5]). The 
state obtained from a given stochastic trajectory X is a 
superposition of 2 R% terms corresponding to all the paths 
along this tree structure where R% denotes the number 
of relaxation events along X (see fig. |9] for an example 
where X has two relaxation jumps, two dephasing jumps 
and two photon jumps). Note that each quantum jump a 
introduces a slow phase factor £ a depending on its type 
and its occurence time (see fig. [SJ. Finally, the state 
\^>%(t)) associated with the quantum trajectory X has 
the form: 



l*l(*)> 



E 

TePath(T) 



& i* m ;(M*))) (B8) 



> 0<t„<t 



where the sum in the r.h.s is performed over all paths T 
going downward the tree structure generated by quantum 
jumps along X (black line in fig. [5]). In this expression, 
m' t = (— l) J I T lm depends on the number J[T] of direc- 
tion changes along the path T. The angle 0r(t) is related 
to the Fresnel angle </>r(i) by <f>r(t) = —m' t 6r{t)/2y/E. 



c. Summing over stochastic trajectories 

With this simple description of the stochastic dynam- 
ics, the evolution of the coherence (IB1I) is of the form: 



n„ 



.(*) 



E < 



{m 1 , ,m'_ ) 



(t) 



(B9) 



where n 



X%}) and \X%! ) states with m[ = ±m e . Explicitely 



( m ' + 



An) 



(t) only involves coherence between the 



n 



{m 1 , ,m'_ ) 
m_i- ,m_ 



(*) 



E 



(p++P-)/2 



p+X Vp+ ] -p- 1 

x (p +' p l ) (f)|A' (p + ) }(X (p r ) | (BIO) 
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fern => S V = 1 f- = -e^/v^ 

(a) (b) (c) 



d. Atomic observables 

Let us assume that the atoms + cavity reduced density 
operator at time t is of the form given by the r.h.s of (|B9|) . 
The Rabi oscillation signal is completely contained in the 
qubit polarization: 

(o*{t)) = e~*J2 J \ (vi?Ht)+T>*?(t)) (B13) 

and clearly involves the evolution of coherences between 
|*+) and |^_). Note however that it does not probe 
the coherence between subspaces associated with differ- 
ent values of p. The qubit coherence is given by: 



FIG. 8: Tree structure associated with a stochastic trajec- 
tory with two photon emission, two pure dephasing quantum 
jumps and two relaxation jumps: (a) Vertices associated with 
the loss of a photon from the cavity (b) Vertices associated 
with a 7r rotation of the qubit along the z§axis (pure dephas- 
ing) (c) Vertices associated with relaxation of the qubit in its 
ground state (relaxation). 




FIG. 9: The diagrammatic representation of a typical quan- 
tum stochastic trajectory 1 showing two photon emissions 
from the cavity, two relaxations and two pure dephasing quan- 
tum jumps. Time would correspond to an horizontal axis 
whereas the vertical axis is the angle of the quasi generalized 
Gea-Banacloche state in the Fresnel plane. A single path T 
along X is singularized (bold line). 



where the coefficient T> , 'Z,/\t) (which will be denoted 

by "D P+ , P _{t) to shorten notation) is given by a path in- 
tegral over all pairs of paths T± downwards all possi- 
ble pairs of quantum trajectories T± of the appropriate 
phase: 

2W(*) - J ■D[T + ,T-]St[T+] x S t [T_r 

x e 2in(y/p+ + l4>T + (t)-y/p- + 14>T_(t)) ^ Bn ) 

Here Stp±] are the products taken along T± of slow 
phases appearing in eqs. (|B5[) and (|B6|) : 



St[T± 



n & 



(B12) 



(t)-V 



(p,p+i) 



(*) 



(p,p+l) 



(*) 



(B14) 



and is sensitive to the coherence between adjacent mul- 
tiplcts. 



Decoherence coefficients 



The main problem is now to compute all decoherence 
coefficients "D p+ P _(t) (which also depends on m± and 
m' ± ). As we shall see, this cannot be done exactly but a 
careful analysis of these sums in the mesoscopic regime 
will show us which trajectories contribute to the Rabi 
oscillation signal. 

Generically, the decoherence coefficient T> p+ , p _ (t) con- 
tains slow phases arising from the phase factors asso- 
ciated with photon emissions and qubit relaxations but 



arising 

r(p±)\ 



also a rapid random phase e 2 " lA "+"- ( ^+- </,r - )(t) 
from difference of the eigenenergies of the states \Xm ± 
Note that in the absence of qubit relaxation and dephas- 
ing [13], this rapid phase is not random and therefore, 
decoherence is completely caused by the slow phases as- 
sociated with photon losses of the cavity. Here, relaxation 
and dephasing quantum jumps of the qubit lead to the 
decay of T> p+tP _(t). Indeed, in the mesoscopic regime, 
averaging of this rapid phase dominates the decay of the 
decoherence coefficient £> p+lP _ (t). To confirm this, let us 
consider the coefficient D p+ , p _ (t) defined by forgetting 
slow phases in 2? p+iP _ (t): 



D p+tP _(t)= / P[T + ,71]e 



2in(y/p+ + l<h+ (t)-y/p-+l&T_ (*)) 

. (B15) 

As we shall see now, this expression can be studied an- 
alytically by adapting renewal theory tec hniq ues already 
used to compute decoherence coefficients [3fJ. 



16 



a. Waiting time probabilities 

The effective statistics of quantum jumps is Poissonian 
(see appendix [5] for relaxation) and statistics of jumps 
of different types are statistically uncorrelated. Equiv- 
alently, the probability distribution for waiting times 
between two jumps of type j is exponential: - 0j( T ) = 
7j e~^' T where 7 cav = fin, j Aeph = 7^/2 and 7 rc i ax = 
7i/2. This exponential factor is taken into account 
through the norm of the states (|A3[) . 



b. Pure dephasing 

In this case, there are no slow phases and no branching. 
No approximation has been made beyond the mesoscopic 
approximation on the dissipationless problem. Hence we 
expect the following results to be valid even in the regime 
of strong dephasing g < 7„. 

The decoherence coefficients 'D^'^ft) are given by 
averages over a telegraphic noise: 

^/°\t) = (J x K x ^ dT ) c/o . 



(B16) 



Here A = g(m+y / p+ + 1 — m-^Jp- + 1) and X(t) = ±1 
is a telegraphic noise characterized by the probability dis- 
tribution of waiting times between two switching iJj(t) = 
(-/ ip /2)e-^ T / 2 and the initial condition X(0) = 1. The 
e (respectively o) symbol specifies that the average is 
performed over noise histories with an even (respectively 
odd) number of switchings between and t (m' ± = m± 
for the e case and m' ± = —m± in the o case). 

Renewal theory provides an elegant analytical solu- 
tion for the decoherence coefficients (|B16[) . Introducing 
Ho(t) = r°° V'( r ) dr and proceeding along the lines of 
[37| , the Laplace transforms of the D^ -* is obtained as: 



L[D[ C) ]( S ) = 



L[Uo](s + iX) 



1 - L[i/>] 


(s - iX)L[ip](s + 


■iX) 








1 - L[i>] 


(a-iA)i[^](a-+ 


iX) 


m(s- 


i\)(l-L[ll>]{8 


+ iX)) 


s + iX 



(B17) 
(B18) 
(B19) 



Performing the inverse Laplace transform gives the ex- 
plicit time dependence. The "odd" coefficient is thus 
given by: 



V {o \t) 



2( S+ -s_) 



(eZ+t-e 3 - 1 ) 



where 



2 



±^JX 2 ~ 1 2 /4. 



(B20) 



(B21) 



The "even" coefficient has a somehow more involved ex- 
pression: 



iX 



(B22) 



two 



Depending on the dimensionlcss coupling g = 2X/j v , 
distinct regimes occur. 

For g > 1, only a single event is enough to spread the 
phase significantly. Decoherence coefficients arc domi- 
nated by noise histories with very low number of switch- 
ing events. For T>£ \ those are histories without any 

switwhing event whereas for T>^\ histories with a sin- 
gle switching event dominate. In this strong dephasing 
limit g > 1, these expressions can be simplified 



x e -iM 



2?f (t) = 



2A 



sin (At) 



(B23) 
(B24) 



In the g> 1 limit, only the "even" contribution survives. 

For g <C 1, a large number of switching events is neces- 
sary to spread the phase A J Q X(t) dr significantly. The 
dephasing time is then much longer than the typical wait- 
ing time. In this regime, we have: 



V 



(e/o) 



(B25) 



leading to a decoherence time of the order of r f v /X 2 ^ 
7~ . Both coefficients are equal since having an even 
or an odd number of switchings makes no difference for 
t>7^- 

In the present context, the coupling constant A de- 
pends on p± and m± . 

For m+ 7^ m_, we focus on the decoherence of a pair 
of generalized Gea-Banacloche states. In this case, for all 
values of p± close to n, the dimensionless coupling is of 
the order of g ~ 2g\/fi,/~f v 3> 1 meaning that we are in 
the ff^l regime. In this case, formulas (|B23[) and (|B24|) 
arc directly relevant and the decoherence rate is given by 

For m + = rri- = m, we focus on the decoherence of a 
single generalized Gea-Banacloche state. In this case, for 
\p + — p_ | < y/E, the dimensionless coupling is then of the 
order of g ~ mg(p + — p^)j^ v \ffi. Its smallest non zero 
value is reached for \p+— p-\ = 1: mg/j^^/E. Therefore, 
for n > n c where n c = (2g/j ip ) 2 ^> 1, the j < 1 limit 
can be reached for small enough \p + —p—\. On the other 
hand, the highest value is g ~ mg/j^ 3> 1. 

As long as n <C n Cl which is the case of practical in- 
terest here, all decoherence coefficients V\(t) associated 
with \p + — p- 1 < \ffi are in the g ^> 1 regime. But for 
n ^> n c , this is no longer the case. In this case, starting 
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with a qubit initially in state \X±), we find that the re- 
duced density operator for the electromagnetic mode is 
of the form: 



i(P++P-)/2 



■e rp +- p - t \p+)(p-\ 



(B26) 

where T p+p _ is a rate bounded from above by 7^ and 
for sufficiently low p+ — p~, given by: 



■ p+, 



The physical interpretation of these results is clear: the 
condition n > n c means that the average number of de- 
phasing jumps within the revival time 4iry/fi,/ g is much 
larger than one. Since each of these jumps is equivalent 
to an echo pulse, the phase of the electromagnetic mode 
has a diffusive motion around zero thus leading to de- 
coherence. In the end, the qubit, incoherent over times 
scales > 7~ 1 due to its coupling with its own environ- 
ment, tends to select specific states of the electromag- 
netic mode. These turn out to be Fock states. The fact 
that the coherence between adjacent Fock states decays 
over a much longer time scale than 7" 1 (see (|B27[) ) re- 
flects the time needed for a small incoherent object (the 
qubit) to decohere the mesoscopic coherent field. 

To summarize, in the regime g <C J v and 1 <C n <C n c , 
the main contribution to the atoms + cavity reduced den- 
sity operator comes from quantum trajectories without 
any quantum jump (decoherence coefficients (t)). 



c. Qubit relaxation 

Being primarily interested in the Rabi oscillation sig- 
nals, we shall focus on the average value (a z (t)). Eq. 
(|BT3)) relates it to the coherences |x| p) )(xi p) | and 



necting to \X_')(X^''\ at time t a 3- valued noise his- 
tory: X{t) £ {0, ±igy/p + 1} so that the pure phase as- 
sociated with the pair of paths the is exp (i J* X(t) dr). 




\X^ ){X_ |. Therefore it is sufficient to compute the \X ( - P ' > )(X 
decoherence coefficient T>^f' p \t). Because of the precise 
form of relaxation quantum jumps, it receives contribu- 
tions from all four coefficients D^f-* (0). As a first guess, 
we might forget about the slow phases and compute ap- 
proximate decoherence coefficients ■ In this case, 
the contribution of each pair of path is a well defined pure 
phase divided by 2 R ^ where i?[T] is the total number of 
relaxation jumps of the underlying quantum trajectory. 

Consider now a specific pair of these paths [7!|_,71]. 
Since we arc dealing with p + = p_, the phase accumu- 
lated between quantum jumps (duration At) is equal 
either to e ±l9v/p+1AT when the two paths 7+ and 71 
are not parallel or to 1 when they are parallel. Thus, 
one may conveniently associate with all [7+,71] con- 



FIG. 10: Plots of all pairs of trajectories involved in the evolu- 
tion of the \X-){X+ \ coherence having exactly two relaxation 
jumps. Relaxation jumps are shown as bifurcations. The + 
and — labels indicate which trajectories are relative to the 
forward and backward Keldysh contours. The extremal tra- 
jectories that dominate (a z (t)) are indicated with thick lines. 



Starting from the coherence (X^ |, there is one 

single pair of paths which connects to the same coher- 
ence and whose associated phase does not depend at all 
on the dates of quantum jumps occuring between t = 
and t. It consists into the pair of extremal paths de- 
picted on fig. [TO! and corresponds to X(t) = g^Jp + 1. 
Any other pairs of paths will lead to a noise X(t) pre- 
senting at least one blip, i.e. a time interval during 
which X(t) 7^ g^/p + 1. Any other pair of paths start- 
ing from \X^)(X^\ or any pair of paths starting from 
the diagonal terms \X { e p) ) (x[ p) | (e 



(?) 



±) and ending on 

will exhibit an associated noise r t— > X{t) 
that is not constant. 

Given a certain noise history r 1— > X(t), we shall now 
sum over all pairs of trajectories [7+, 71] associated with 
a given function r 1— > X(t). One of them has a min- 
imal number of relaxation jumps precisely occuring at 
the dates where X{t) jumps. But of all histories having 
more relaxation jumps occuring between these specific 
dates must also be taken into account (this is an impor- 
tant difference with the case of pure dephasing studied in 
the previous paragraph). As shown in appendix |A1 the 
relaxation jumps stochastic process can be approximated 
by a renewal process with exponential waiting time dis- 
tribution: iPr(t) = 7e" 7T where 7 = 71/2. Taking into 
account the factor 1/2 associated with each relaxation 
jumps, the resulting integration measure for intermedi- 
ates times (t±, . . . , t p ) where < ii < . . . < t p < t is 
given by dfi(ti, ...,t p ) = e" 7 '/ 2 (j/2)P dt ± . . . dt p . Note 
that the exponential factor e -7 * present in the waiting 
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time distribution is partially compensated by the sum- 
mation of all products of 1/2 factors associated with re- 
laxation jumps occuring at times different from the tjs. 
The 1/2 factors associated with jumps occuring at times 
tj are taken into account through ( r y/2) p . 

This counting argument makes it clear that the con- 
tribution of paths associated with a non constant noise 
history will, comparatively to the case X(t) = X(t), in- 
volve integrals of phases of the form e l9T%/p+1 or e 2l9T ^ p+1 
over t with measure /2)dr, thus leading to a factor pro- 
portional to 7/<7y / p + 1 which, in the mesoscopic regime, 
is or the order of j/g^/E. Besides this, we are in the 
strong coupling regime and therefore 7/1? V™ is always 
much smaller than one. Thus, these contributions vanish 
in the regime considered in this paper. 

Finally, our analysis shows that the dominant contri- 
bution to Pjf^(t) in (jBlip comes from trajectories T± 
for which the Fresnel angle (f>x ± is monotonous in time. 
Of course, our argumentation ignored the slow phases 
that appear in (|B12[) but it can be checked by explicit 
computation that our conclusion remains valid for pairs 
of paths with one or two jumps. 

Retaining only the contribution of the pair of extremal 
paths in (|B11[) and taking slow phases into account, 
the argumentation of [l(| can then be straightforwardly 
adapted to resum their contribution along these trajec- 
tories. This immcdiatly leads to (f2"B")) where the 1/2 pref- 
actor in front of e^ A< ^ T ' comes from the 1/2 factor as- 
sociated with each relaxation jump. 



As a final comment, the above derivation also makes 
it clear that depending on the physical quantity we are 
interested in, other pairs of trajectories will contribute. 
As an example, computing the probability for having the 
qubit in the |+) or the |— ) involves coefficients T>^^ 

V { ™ ] which are dominated by the sum over pair of paths 
such that T + — 71 . 



d. Summing over all dissipative processes 

In the presence of dephasing, the extremal trajectories 
are still the dominant contribution to (a z (t)) and the ef- 
fect of dephasing is an extra e -7 *''/ 2 factor. Note that 
these corresponds to quantum histories without any de- 
phasing quantum jump. These are expected to provide 
the dominant contribution in the low 7^ limit. 

It follows from the previous analysis that the main con- 
tribution to (<7 z (t)) comes from pairs of extremal paths. 
Since photon losses are independant processes from relax- 
ation jumps, their contribution can be computed using 
the formalism developped in Ref. [l(| . Obviously, it fac- 
tors in front of the relaxation contribution. For Kt <C 1, 
this enables to treat the regime where photon losses lead- 
ing to decohcrence of the generalized Gea-Banaclochc 
states occur at a much higher rate than relaxation and 
pure dephasing nn 3> 7i, v - 
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